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Abstract. We present a Monte Carlo study of dijet angular distributions at y/s — 14 TeV. First we perform 
a next-to-leading order QCD study; we calculate the distributions in four different bins of dijet invariant 
mass using different Monte Carlo programs and different jet algorithms, and we also investigate the sys- 
tematic uncertainties coming from the choice of the parton distribution functions and the renormalization 
and factorization scales. In the second part of this paper, we present the effects on the distributions com- 
ing from a model including gravitational scattering and black hole formation in a world with large extra 
dimensions. Assuming a 25% systematic uncertainty, we report a discovery potential for the mass bin 
1 < Mjj < 2 TeV at 10 pb"^ integrated luminosity. 
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1 Introduction 

Jet produetion is the most dominant hard proeess in hadron 
collision experiments. While jets are background for many 
new physics searches, jets can also be used as a signal, 
both for probing QCD and for physics beyond the Stan- 
dard Model. Because of their rich abundance, many jet 
studies can be performed with little integrated luminos- 
ity. 

The dijet angular distribution between the two hardest 
jets in the event has shown to be a very useful measure- 
ment [UElElll] ; at low integrated luminosity it is a good 
tool to probe QCD, while with more statistics, a search 
for new physics, such as effects coming from large extra 
dimensions, becomes possible. 

This paper is about dijet angular distributions at ^/i = 
14 TeV. First we will perform a QCD study; we will cal- 
culate the distributions up to NLO and make an estimate 
of the systematic uncertainties. In the second part we 
will study the distributions in a scenario with large extra 
dimensions, including gravitational scattering and black 
hole formation. 



2 Dijet angular distributions 

In proton-proton collisions, QCD will mainly manifest it- 
self in t-channel gluon exchanges, giving rise to a cross 
section which is approximately inversely proportional to 



with i 



da 

— - cx 
dt 



£2 



(1) 



;(1 — cos 0), Q being the scattering angle in the 



center of mass frame of the colliding partons. This corre- 
sponds to Rutherford scattering for constant s, i.e. a scat- 
tering process that peaks at small angles: da j dico^ ff) cx 



sm 



Experimentally, the study of the angular behavior is 
done using the variable x, which is defined as x = expdjyi — 
rii\) = exp(2|r/*|), with r\x and r\2 the pseudorapidities of 
the two hardest jets in the event and rf = ^(jyi — 112)- At 
lowest order, i.e. for a 2 —> 2 process, the pseudorapid- 
ity of the two jets (partons) in the center of mass frame 
is given by ±77* and can be related to the scattering an- 
gle 9 = arccos(tanh7y*), and therefore also the following 
expression for x as a function of 6 holds: 



1 



I cos 6*1 
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Using the approximation that x cx s/i and keeping s fixed, 
the cross section in equation ([1]) can now be rewritten as a 
function of x and turns out to be approximately constant: 



— cx — (s fixed) 
dx s 



(3) 



To get the result at the hadron level, the above cross sec- 
tion needs to be multiplied with the parton distribution 
functions and integrated over the momentum fractions: 



da 



dxi / dx2fl{xi,Q'^)f2{x2,Q" 



I? 



(4) 



Keeping s fixed implies that to a good approximation, 
also the product of the two parton distribution functions 
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is fixed (up to logarithmic scaling variations with Q^), so 
that after integration, the cross section da/dx becomes 
constant as well. This means that, binned in the dijet in- 
variant mass Mjj = Mjj at LO), the angular dis- 
tribution da/dx versus x is approximately flat for QCD. 
On the other hand, new physics is very often an s-channel 
resonance or a process characterized by a certain mass 
threshold. These processes are usually more isotropic than 
QCD, with a cross section that tends to be flat in cos 9, 
which causes the dijet angular distributions to peak at low 
X- Therefore, in order to gain most knowledge from dijet 
angular distributions, it is necessary to study them in bins 
ofM,,. 

The following four mass bins were chosen for y^=14 
TeV: 0.5 < Af^j < 1 TeV, 1 < Mjj < 2 TeV, 2 < Mjj < 3 
TeV and 3 TeV < Mjj. 

In experimental conditions the measurable phase space 
will limit the maximum x value one can reach. If r/max is 
the maximum pseudorapidity for which jets are still fully 
seen by the detector, then the following two orthogonal 
selection cuts need to be made in order to measure the an- 
gular distribution without acceptance losses coming from 
the limited pseudorapidity range. 



1 771 - 772 I < 277,nax - C 

\m + m\ < c 



(5) 
(6) 



For example the ATLAS [5] calorimeters can measure jets 
fully up to ?7inax ~ 4. The choice of c is made by consider- 
ing the trade-off between having enough statistics passing 
selection cut ^ and the measurable x range, and for our 
study c = 1.5 has turned out to be a reasonable compro- 
mise. With this choice for the value of c and with Tymax = 4, 
the angular distributions can be measured up to Xmax ~ 
600. In case we are only interested in measuring up to 
Xmax ~ 100 (e.g. for new physics searches), we can limit 
ourselves to Ty^ax = 3.1. 

At lowest order, the relation between Mjj, x a-^d the 
transverse momentum pT is the following: 



(7) 



Table 1. Values of for 4 different mass bins and 2 values 

of Xmax- 



J-]j,nm\ < Mjj 



< M. 



jj,max 



The selection cuts on dijet mass (A/- 

and X (x < Xmax), will determine the minimum the two 
highest jets in the event need to have in order to pass the 
selection cuts: 



Mi 



PT,min,LO 



(8) 



NLO contributions will lower the minimum transverse mo- 
mentum with a factor a/S, as a consequence of the fact that 
the second highest jet in the event can never have a trans- 
verse momentum less than half the transverse momentum 
of the highest jet. So. up to NLO: 



M. 



PT.min 



V^(\/ Xmax ^~ l/-\/Xmax) 



(9) 



Equation ^ will be used as a pT-selection cut. Table [T] 
summarizes this for the different mass bins and for two 
different values of Xmax ■ 



Mass bin 


Xmax = 100 


Xmax = 600 




pT^i„ (GeV) 


pT„i„ (GeV) 


0.5 < Mjj < 1 TeV 


35 


14 


1 < Mjj < 2 TeV 


70 


28 


2 < A4jj < 3 TeV 


140 


57 


3 TeV < Mjj 


210 


86 



3 QCD calculations 

Two programs are available for NLO jet calculations: JE- 
TRAD [g and NLOJET-h-h [7]. The programs use a con- 
ceptually very different approach; JETRAD uses the phase 
space slicing technique, while NLOJET-I--I- applies the 
Catani-Seymour dipole subtraction scheme [5] with some 
modifications introduced because of computational rea- 
sons. 

While the considerations in the previous section are 
at the parton level, experiments have to deal with jets 
made of hadrons, which are the result of three major 
steps: hard interaction, parton showering and hadroniza- 
tion. Jets carry the memory of the hard interaction and 
a good jet-finding algorithm can exhibit this information 
in an infrared and coUinear safe way. For many years, 
the longitudinally-invariant /cT-clustering algorithm has 
shown to be infrared and coUinear safe [9lll0j. and it will 
be used in this study. Recently, SISCone [11], a practical 
infrared safe cone algorithm, has become available and 
we will also show results using this algorithm. Although 
seeded cone algorithms are not infrared stable, we will also 
present results with a seeded cone algorithm because they 
are still frequently used in experiments. More precisely, 
we will use a seeded, iterative cone with progressive re- 
moval that comes with JETRAD and was brought into 
use by the CERN UAl collaboration [12]. NLOJET-h-F 
comes with an exact seedless cone algorithm which was 
first proposed in [13j and is finding all stable cones in a 
given configuration, ensuring infrared safety. We will also 
show results using this algorithm. 

Note that the final state for a NLO order calculation 
at the matrix element level contains at most three partons 
and that we consider only those events with at least two 
jets. Therefore most problems emerging from infrared and 
coUinear instabilities are absent in our calculations. The 
major difference between the previously mentioned cone 
algorithms is that the JETRAD seeded cone clusters two 
nearby partons if their separation in (77, (p) is less than the 
cone radius R, while the NLOJET-|--|- seedless cone and 
SISCone only do so if the cone containing both partons 
is stable, i.e. if their separation is smaller than R{1 + z), 
with z — pt,2/pt,i and pt,2 < Pt,i [Hj- 

Figure [1] shows the angular distribution for the mass 
bin 1 < Mjj < 2 TeV, calculated both with JETRAD and 
NLOJET-f +, using an inclusive kx algorithm with radius 
parameter R = 1.0. JETRAD uses a different parametriza- 
tion of the strong coupling constant than NL0JET+4-, 
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Fig. 1. Comparing JETRAD with NLOJET++ for the mass bin 1 < Aljj < 2 TeV. The calculations are done at NLO using 
an inclusive kr algorithm with R = 1.0. 



which explains the difference between the curves to a large 
extent. To illustrate this, we show in figure [5] the angular 
distributions with as kept constant at 0.1. 

NLOJET++ (coded in C++) has the advantage over 
JETRAD (coded in fortran) that it can be combined with 
more modern jet algorithms, such as SISCone and an exact 
seedless cone. But it has the disadvantage that the angular 
distributions have a statistical error that is not homoge- 
neous over the whole x range, but is instead increasing 
with X, which is not the case for a JETRAD calculation, 
as can be seen from figure [1] 

Figure [3] compares calculations at the Born (lowest or- 
der) level with next- to- leading order calculations, done 
with JETRAD, for the four different mass bins and for 
X < 600. ALL NLO distributions have been calculated 
with the next-to-leading order CTEQ6M parton distribu- 
tion function (PDF) [14] , while the leading order distribu- 
tions have been calculated both with the CTEQ6M PDF 
and with the leading order CTEQ6L1 PDF. Two different 
jet algorithms are used for the NLO calculations; the JE- 
TRAD seeded cone algorithm with radius R=0.7, and an 
inclusive kr algorithm with R = 1.0. Note that in a LO 
parton level calculation, the outgoing partons are back- 
to-back, so that a jet algorithm is redundant. The NLO 
angular distributions with the two different jet algorithms 
tend to have the same shape, but differ in absolute nor- 
malization. 



From this figure, we also note that the NLO calcula- 
tions are flatter than the Born calculations, especially at 
high X values, which can be explained to a large extent by 
the fact that the running of as with pT (or equivalently 
x) has more effect on a LO than a NLO calculation. To 
illustrate this, we have plotted in figure [4] four different 
LO calculations in the mass bin 1 < Mjj < 2 TeV; the 
dashed blue curve is the same Born calculation as the one 
presented in the top right plot of figure [3l For this calcu- 
lation, the value of as and of the factorization scale (fJ-p) 
have been varied according to the pT of the hardest jet. 
The other three curves show what happens if as and/or 
fiF have been kept constant. The grey dotted curve is a 
calculation with as = 0.1 and /i_F = 100 GeV and is much 
flatter compared to the blue one. The red curve has been 
calculated with as = 0.1 and no running oi ^p, and the 
green one with = 100 GcV and no running of as- There 
is a much bigger effect on the distributions from keeping 
as flxed than from ^p- 

But these are observations at the Born level. At NLO, 
the sensitivity due to scale variations is reduced as can be 
observed in figure [5l This is because perturbation theory 
tells us that an all-order calculation should not depend on 
the renormalization scale at all, and therefore, compared 
to the LO calculation, a NLO calculation is more stable 
against the running of as- Furthermore, a NLO order cal- 
culation is using a NLO expansion of as , and the running 
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Dijet angular distributions 1 < Mjj < 2 TeV, ag = 0.1 
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Fig. 2. Comparing JETRAD witii NLOJET++ for tlie mass bin 1 < Mjj < 2 TeV and for as = 0.1. The calculations are done 
at NLO using an inclusive kr algorithm with R — 1.0. 



of as at NLO is less pronounced than at LO. Figures [Hand 
[5] clearly illustrate the need for a NLO order calculation 
with the running of and enabled. 

This is also confirmed by calculations done with NLO- 
JET++. In figure [5] we present the angular distributions, 
done with NLOJET++, for the same mass bins, using the 
NLOJET++ seedless cone algorithm with R = 0.7 and 
overlap 0.5, the SISCone algorithm with R=G.7 and over- 
lap 0.75, and an inclusive kr algorithm with R = 0.6. We 
show the distributions for x < 100 only and, compared to 
the big variations between the LO and NLO calculations 
observed at high x (x > 100) in figure [3l the difference 
between LO and NLO is much less at small values of x- 

In figure [7] we show the ratio of the NLO cross section 
calculated with CTEQ6M over the LO one calculated with 
CTEQ6L1, for x < 100 and for different jet algorithms, 
a quantity which is often called the k-factor in literature. 
The k-factor varies around 1 with a deviation of 30% at 
the most. 

The calculations above have been done with the CTEQ6M 
and CTEQ6L1 parton distribution functions and with a 
normalization and factorization scale chosen to be the 
transverse momentum of the hardest jet. Uncertainties 
coming from parton distribution functions (PDFs) and 
the choice of renormalization {(ir) and factorization (fip) 
scale will contribute to a systematic error. We will investi- 
gate them using JETRAD with an inclusive kr algorithm 



with R = 1.0. The exact scale and PDF uncertainties may 
vary with different jet algorithms and cone sizes. 

Figure [8] shows the angular distributions for the mass 
bin 1 < Aljj < 2 TeV for three different PDF-sets, namely 
CTEQ6M, CTEQ66 ^ and MSTW2008NLO [B]. The 
distributions differ mainly in absolute normalization, and 
less in shape. Normalizing the distributions to unit area 
make them differ no more than 3% over the whole x range, 
as can be seen in figure [9l 

To further examine the uncertainties coming from PDFs, 
we have calculated the angular distributions for all 44 er- 
ror members of the CTEQ66 PDF and applied the Master 
Equation suggested in [T7] to deduct a positive and nega- 
tive uncertainty on a quantity X: 



AX+ = 

max 



\ 



N 

E 



[max(X+ - Xo, Xr - Xo, OW (10) 



AX+„ = 



\ 



N 



J2[max{Xo - X+,Xo - Xr,0)]^ (U) 



AX^ adds in quadrature the PDF error contributions that 
lead to an increase in the observable X, and AX~ the 
PDF error contributions that lead to a decrease. Using 
this formula, the error on the dijet angular distribution 
calculated with the central PDF member, in the mass bin 
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Fig. 3. LO and NLO angular distributions calculated with JETRAD for 4 different mass bins. The NLO calculations are done 
with two different jet algorithms: the JETRAD seeded cone algorithm with R=0.7, and an inclusive fcr algorithm with R = 
1.0. All NLO curves are calculated with the CTEQ6M PDF, while the LO distributions are calculated with both the CTEQ6M 
and the CTEQ6L1 PDF. In the bottom two plots, the LO CTEQ6L1 and LO CTEQ6M curves are overlapping. 
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Fig. 5. Illustrating the influence of the running of Os and the variation of /if on NLO distributions. The calculations are done 
with JETRAD and an inclusive algorithm with R = 1.0. 
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Fig. 6. LO and NLO angular distributions calculated with NLOJET++ for 4 different mass bins. The NLO calculations are 
done with three different jet algorithms: an inclusive kx algorithm with R = 0.6, SISCone with R = 0.7, and the NLOJET++ 
seedless cone algorithm with R — 0.7. 
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Fig. 7. Ratio of aNLo(CTEQ6M) over aLo(CTEQ6Ll) for 4 different mass bins calculated with NLOJET++. 




Fig. 8. Angular distributions with different PDF-sets for the mass bin 1 < Mjj < 2 TeV. The calculations are done with 
JETRAD and an inclusive kr algorithm with R = 1.0. 
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Fig. 9. Difference caused by tlie use of different PDF-sets on botli normalized and unnormalized dijet angular distributions for 
the mass bin 1 < Mjj < 2 TeV. 
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Fig. 10. Central member of the CTEQ66 PDF, together with its uncertainty band for the mass bin 1 < Mjj < 2 TeV. The 
calculations are done with JETRAD and an inclusive kr algorithm with R = 1.0. 



1 < Mjj < 2 TeV is plotted in figure [10] as a blue band 
around the central member. The uncertainties given in 
terms of percentage are shown in figure [TT] 

The choice of fin and fip will also influence the dis- 
tributions. We have studied this by letting fi^ and fip 
vary independently between 0.5, 1 and 2 times the trans- 
verse momentum of the hardest jet, resulting in 9 different 
distributions in total. Figures [T2l and [T3l summarize the re- 
sults for the mass bin 1 < Mjj < 2 TcV. The other mass 
bins have similar results. For all these figures, the plots 
at the left show the distributions for x < lOOi while the 



plots at the right go up to x = 600. The variables r and / 
used in the figures, denote the fraction of the transverse 
momentum of the highest jet at which respectively the 
renormalization and factorization scale are evaluated. In 
figure \Vl[ fip is kept constant, while fin is varied, which 
infiuences both the shape and the normalization of the 
distributions. The effect is small at low x but increases 
drastically with increasing x values. The plots in figure [T51 
all have fixed (jLr and varying fip, which causes a change 
rather in normalization and not so much in shape. 
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Fig. 11. Uncertainty on the angular distribution calculated with the central member of the CTEQ66 PDF, in the mass bin 
1 < Mjj < 2 TeV coming from the intrinsic uncertainty of the PDF. The + error leads to an increase of the cross section, 
while the — error leads to a decrease (see equation in the text). The uncertainties are shown on the curves both with and 
without normalization to unit area 1 < x < 100. 




Fig. 12. Influence of the choice of fiR for the mass bin 1 < Mjj < 2 TeV. The variables r and / denote the fraction of the 
transverse momentum of the highest jet at which respectively the renormalization and factorization scale are evaluated. Left: 
X < 100, right; x < 600. The calculations are done with JETRAD and an inclusive kr algorithm with R = 1.0. 



The difference between these distributions is an esti- 
mate of the uncertainty coming from the choice of ^xr and 
/j,p. In figure [TH we show the uncertainty on the distribu- 
tions both with and without normahzation to unit area 
1 < X < 100. NormaHzing the distributions reduces the 
error drastically; averaged over the whole x range, the 
uncertainty on the normalized distribution is 9%, and is 
never exceeding 20%. 



Figure [15] shows the combination in quadrature of the 
uncertainties coming from the choice of renormalization 
and factorization scale (figure [T4|l. together with the in- 
trinsic uncertainty from the CTEQ66 PDF (figure [TT|) . 
drawn as an error band around the calculation done with 
nominal values (central CTEQ66 member and fxu = = 
Pt of the highest jet). Both the distributions with (plot 
right) and without (plot left) normalization to unit area 
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Fig. 13. Influence of the choice of hf for the mass bin 1 < Mjj < 2 TeV for different values of hr. The variables r and / 
denote the fraction of the transverse momentum of the highest jet at which respectively the renormalization and factorization 
scale are evaluated. Left: x < 100, right: x < 600. The calculations are done with JETRAD and an inclusive /ct algorithm with 
R = 1.0. 



1 < X < 100 E^re shown. In both plots, the error band 
is dominated by the error coming from the choice of the 
scales, more precisely, the renormalization scale introduces 
the major uncertainty. When normalized to unit area 1 < 
X < 100, the combined error docs not exceed 20% over 
the whole x range. 

Above calculations have been done at the parton level 
of the hard interaction, without showering, multiple in- 
teractions or hadronization. In figure [11] wc show what 



happens if wc turn on these processes. In the top fig- 
ure, wc have plotted angular distributions calculated with 
PYTHIA 6.410 [IZj. The dashed gray curve is a calcula- 
tion at the Born level, without switching on any of the pro- 
cesses that were previously mentioned. Initial state radia- 
tion (ISR) and primordial kx (prim kT) have been turned 
on for the dashed red curve and the full blue curve is a 
calculation with final state radiation (FSR) included (but 
initial state radiation and primordial kx turned off). As 
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Fig. 14. Systematic uncertainty coming from the choice of jJ.R,F for the mass bin 1 < Mjj < 2 TeV, both on normahzed and 
not normalized distributions. 
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Fig. 15. Black line: calculation done with nominal central CTEQ66 member and fiR = mup = Pt of the highest jet, both 
with (plot right) and without (plot left) normalization to unit area 1 < x < 100. Blue band: error band from combining the 
uncertainties coming from the choice of renormalization and factorization scale (figure[l4}, together with the intrinsic uncertainty 
from the CTEQ66 PDF (figure [TTJ- 



can be observed in the figure, final state radiation causes 
losses out of the jet cone. The dashed magenta curve cov- 
ers both initial and final state showers, primordial fcy and 
multiple interactions. Finally, the full green curve includes 
all these processes and hadronization. The calculations are 
done with the PYCELL cluster routine from the PYTHIA 
library and the PYTHIA default settings for initial and 
final state radiation, multiple interactions and hadroniza- 
tion. For this particular choice of jet algorithm and set- 
tings, the difference between the calculation at the Born 
level and the one at the hadron level is rather small, as can 
be seen from their ratio, shown in the bottom plot of figure 



1161 The corrections may vary for different jet algorithms 
and PYTHIA settings. 

4 Gravitational scattering and black hole 
formation in large extra dimensions 

An attempt to resolve the hierarchy problem in particle 
physics, i.e. why the electroweak scale is so many orders 
of magnitude lower than the observed Planck scale, was 
done by Arkani-Hamed, Dimopoulos and Dvali in the late 
nineties P^SIITO] . Their so-called ADD model assumes the 
existence of large extra spatial dimensions in which gravity 
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Fig. 16. Including primordial kr, initial and final state radiation, multiple interactions and hadronization in a LO calculation. 
The calculations are done using PYTHIA and the PYCELL cluster routine. 



is allowed to propagate, while the SM fields are confined 
to a four-dimensional membrane. The fact that the extra 
dimensions have not been observed yet can be explained 
by assuming that they are compactified with a (common) 
compactification radius R. This causes the fundamental 
Planck scale Mp to be much smaller than the observed 4- 
dimensional one, Mp4, i.e. for n extra dimensions the re- 
lation Mp^ ~ Mp~^^R" holds. Experiments conducted to 
test the gravitational inverse-square law have found that 
an extra dimension must have a size i? < 44 pim . 

For a fundamental Planck scale of around 1 TeV, the 
ADD model predicts the production of extra dimensional 
black holes at the LHC. Besides black holes, also processes 
involving the exchange of virtual Kaluza-Klein (KK) modes, 
with gravitational scattering of hard partons as dominant 
process, will be present. 

Effects from Black Hole production and gravitational 
scattering on standard jet observables such as the px- 



spectrum and the dijet invariant mass spectrum were in- 
vestigated in [3T]. The study was based on an effective 
field description of gravitational scattering done in 
[23l[24] . This section will continue this study by elaborat- 
ing on dijet angular distributions. First we will give a brief 
summary of the theory, then we will show results. 

4.1 Gravitational scattering 

In a gravitational scattering event, quantized KK modes 
will occur as intermediate states and, in order to calculate 
the scattering amplitude, the sum over all modes has to 
be made. This sum (which can be turned into an integral) 
diverges, but can be rendered finite by introducing a phys- 
ical cut-off of the KK tower to be summed over. A cut-off 
implied from a narrow (small compared to the compactifi- 
cation radius R) width of the four-dimensional membrane, 
assuming a Gaussian extension of the standard model field 
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densities into the bulk, was proposed in |24| and a calcu- 
lation of the i-channel amplitude was made. The finite 
width suppresses the coupling to high-mass KK modes 



with a factor e ™ Z*^" 



being the mass of the mode 



and 1/Ms the width of the four-dimensional membrane, 
and an effective propagator can be derived. The scatter- 
ing process now depends on three energy scales, namely 
the collision energy Vs, the fundamental Planck scale Mp 
and the inverse of the membrane width M^. 

For ^ Ms, forward scattering via the t-channel 
dominates and an all-order eikonal calculation was done 
to ensure unitarity. Contributions from multi-loop ladder 
diagrams exponentiate and the all order eikonal amplitude 
is given by 



1), 



(12) 



with b the impact parameter and x(&) the eikonal scatter- 
ing phase: 



X(6) = 1 



(2^ 



ABorn(-fci), (13) 



The Born term used in equation is dominated by 
small values of (t/s), i.e. small angle scattering, and is 
— neglecting spins of colliding partons — given by: 



A 



Born 



/ 7\i s \ " 



\MpJ Ml-t 



(14) 



On the other hand, when the exchanged momentum 
is small compared to Ms, i.e. for ^f^ < Ms, s-,t- and 
u-channels are equally efficient and scattering is almost 
isotropic. KK propagators in the f-channel can be replaced 
by vertex factors so that the exchange of KK modes cor- 
responds to a contact interaction. A geometric series of 
ladder diagrams is obtained and summed over so that uni- 
tarity is constrained: 



A 



with X 



ladders 



(In + in) . The Born term — neglecting 
spins of the colliding partons — can be written in terms of 
an effective Planck mass Mes- 



(15) 



probability of a collision between two partons with iden- 
tical flavor, spin and color. 

Note that gravitational scattering will become weaker 
for smaller values of Ms (i.e. wider membranes), because 
the KK-modes with a mass above Ms are suppressed. 



4.2 Black holes 

The Schwarzschild radius, introduced by Schwarzschild for 
static non-spinning massive objects, can be generalized to 
4-|-n dimensions 1251 : 



1 ,Mbh- 
Mp ^ Mp ' 



(18) 



In hadron collisions, black holes can be formed when the 
interacting partons come closer than twice . This means 
that in the rest frame of the incoming partons their longi- 
tudinal wavelength A/ ex 2/ \/I and transverse wavelength 
A_L oc 1/pt need to be smaller than Vg. This implies a 
minimum on the black hole mass, which lies in the neigh- 
borhood of the Planck scale [26] : 



M, 



min,l 



= A/p( 



(20?) 



n+l. 



2) 



)^ 



(19) 



Another limit on the black hole's mass comes from 
the existence of the finite width of the membrane, be- 
cause, in the approximation of a narrow width, the mem- 
brane cannot be more extended than the black hole itself. 
This means that the Schwarzschild radius should not 
be smaller than 1/Ms, which is only true for masses above 

Afmin,2: 



M, 



min.2 



8^(3±Ii)M^^ 



(20) 



The maximum of (|20p and (jl9p is the minimum black 
hole mass possibly created. For n=4 extra dimensions and 
for Ms/Mp > 1.12 {Ms/Mp > 1.05 for n=6), the min- 
imum black hole mass equals (|19p and depends on Mp 
only. In the complementary region the minimum black 
hole mass is described by and goes as Mp/{Ms/Mp)^ 
{Mp/{Ms/Mpy for n=6). In the latter case, small values 
of Ms, i.e. larger values of the membrane width, will pre- 
vent black holes from being created. 



Ae 



with: 



M,. 



1 (n-2)2"7r— 
2^ M"-2 



(16) 4.3 Results 



(17) 



The s-channel ladder diagrams (rotated by 7r/2) are uni- 
tarized in a similar way. KK modes can also be produced 
on shell in this channel, but with a lifetime too long to 
detect. 

The relevance of the u-channel contribution is sup- 
pressed for proton-proton collisions because of the low 



Gravitational scattering was implemented in the PYTHIA 
6.410 event generator and combined with the CHARYB- 
DIS [55] black hole generator, in a similar way as in |21| . 
The calculations have been done using an inclusive kr al- 
gorithm with separation parameter 1.0. We have used the 
k-factor which was derived in the previous section (see 
figure [7]) to scale the PYTHIA QCD distributions up to 
NLO. 

As discussed in the previous sections, the model pa- 
rameters determining the phenomenology are Afp, n and 
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Table 2. Different parameter sets 



Name 


Meff (TeV) 


n 


Ms/Mp 


Mp (TeV) 


Ms (TeV) 


CI 


1.0 


6 


1.0 


0.282 


0.282 


C2 


1.0 


6 


2.0 


0.564 


1.128 


C3 


1.0 


6 


4.0 


1.128 


4.513 


C4 


0.5 


6 


8.0 


1.128 


9.027 


C5 


1.0 


6 


8.0 


2.257 


18.05 


C6 


1.0 


4 


4.0 


1.263 


5.053 



Mg. An equivalent set of parameters is M^s (equation 
(fT7|) ). Ms/Mp and n, and we will work with the latter 
one. Table [2] shows a few combinations. All parameter 
choices in tabic [2] have been tested at a Tcvatron energy 
of Vs = 1.8 TeV for the mass bin 635 GeV < Mjj and 
for X ^ 20, and turned out to be consistent with the dijet 
angular distribution measurements done by the CDF [3] 
and DO collaboration [4|. 

To get a feeling of the impact of the parameters, we 
have performed a few runs with PYTHIA with parton 
showers, multiple interactions and hadronization turned 
off. Figures [T71 to show the dijet angular distributions, 
with and without normalizing to unit area 1 < x < 100, 
for some of the parameter sets defined in table [2l Pre- 
vious experiments have shown that normalizing the dis- 
tributions reduces the experimental error 00. Further- 
more, in the previous section we have demonstrated that 
the theoretical error on normalized distributions is smaller 
as well (figure [T5)) . We have plotted the different contribu- 
tions — i.e. gravitational scattering (GS), black holes (BH) 
and QCD — separately, as well as their sum. As with most 
new physics, the biggest effects from gravitational scatter- 
ing and black hole formation show up at low x (high pr) 
values. How much each process contributes to the total 
cross section, depends on the parameter settings and the 
mass bin. We will make a few quantitative observations. 

Increasing M^jMp while keeping M^^ fixed, will result 
in an increase of gravitational scattering and a decrease of 
black hole formation. The effect is clearly visible when set 
CI (figure [T7)) is compared with set C3 (figure [T5|) for the 
mass bin 1 < Mjj < 2 TeV. For CI, the total cross section 
at low X values, is dominated by black holes, while grav- 
itational scattering is of no importance. But this changes 
drastically for C3; black holes have almost completely dis- 
appeared for C3, but the gravitational scattering cross sec- 
tion has increased by two orders of magnitude. QCD is still 
dominating in this mass bin, but gravitational scattering 
starts to dominate the cross section for mass bins > 2 
TeV, as can be seen in figure fT9l 

Compared to C3, C4 has the same fundamental Planck 
scale and number of extra dimensions (1 TcV and 6, re- 
spectively) , but a different width of the membrane (9 TeV 
vs 4.5 TeV), and this causes effects from gravitational scat- 
tering to set in at lower mass values; from figure [20l it is 
observed that gravitational scattering dominates the cross 
section from 1 TeV onwards, which is not the case for C3 
(figure [Til). 



Figure [21] examines what happens if we double the 
effective and fundamental Planck scale; we have plotted 
the angular distribution in the mass bin 3 TeV < Mjj 
TeV for C4 (figure [21] left) and C5 (figure [21] right). As 
expected, gravitational scattering becomes weaker for C5, 
but it is still strongly dominating QCD. The cross section 
for black hole formation has decreased much more than 
the one for gravitational scattering. 

Finally also the number of dimensions matters. The 
model predicts that the more extra dimensions there are 
the bigger the effects are. Figure gives the angular dis- 
tributions in the mass bin 3 TeV < Mjj TeV for 6 (param- 
eter set C3, plot at the left) and 4 (parameter set C6, plot 
at the right) extra dimensions. The difference between 6 
and 4 extra dimensions is very small. 

Note that when gravitational scattering is the most 
important contribution, deviations from the QCD cross 
section are still visible at large x values, which is not the 
case for black holes. But because at large x values gravita- 
tional scattering is mainly a t-channel process, the shape 
of the graviational scattering curve is close to the QCD 
one, and the difference between QCD and gravitational 
scattering disappears in the normalized distributions. 

Deviations from QCD arc most visible at low x values, 
and a summary of the above observations is done in table 
[3] here we give for a selection of parameter settings and 
mass bins, the integrated cross sections between x = 1 
X = 50 for QCD (ctqcd), gravitational scattering (ctgs) 
and black holes (ctbh) separately. 

We will now compare the different new physics scenar- 
ios with QCD for a given integrated luminosity, so that 
we can establish a discovery potential. We use the dis- 
tributions that are normalized to % < 50 and perform a 
chi-square (x^) test between them. 

We have used the following recipe. For a given inte- 
grated luminosity, we consider the normalized distribu- 
tions {dNQci)/dx)/NQCB and ((iA^totai/rfx)/^totai, with 
Nqcd and iVtotai respectively the number of QCD and 
total (= QCD -I- GS -I- BH) events. We perform a chi- 
square (x^) test between these distributions to test the 
null hypothesis that (c?A^totai/rfx)/A^totai follows the QCD 
distribution. We use both a statistical and systematic un- 
certainty for the calculation of x^- 

X^= E [ I (21) 

all bins I 

where Nqcba cind A^totai,i ^rc the number of QCD and to- 
tal events respectively in bin i. The statistical error 5stat,i 
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Fig. 17. Angular distributions for CI (see table [2)| for the mass bin 1 < Mjj < 2 TeV. Left: cross section in mb, right: cross 
section normalized to unit area 1 < x < 100. 
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Fig. 18. Angular distributions for C3 (see table [2| for the mass bin 1 < Mjj < 2 TeV. Left: cross section in mb, right: cross 
section normalized to unit area 1 < x < 100. 



Table 3. Different parameter sets and the relevance of GS and BH in several mass bins. 



Name 


Meff (TeV) 


n 


Ms/Mp 


Mp 


Ms 


Mass bin 


CQCU 




0-GS (nb) 


0"BH 












(TeV) 


(TeV) 


(TeV) 


(nb) 




(nb) 


(nb) 




CI 


1.0 


6 


1.0 


0.282 


0.282 


[1,2] 


7.23 




8.23 10"^ 


113 




C2 


1.0 


6 


2.0 


0.564 


1.128 


[1,2] 


7.23 




4.74 10"^ 


20.5 




C3 


1.0 


6 


4.0 


1.128 


4.513 


[1,2] 


7.23 




7.01 10-1 


1.01 10" 


1 


C3 


1.0 


6 


4.0 


1.128 


4.513 


[3, 14] 


2.30 10" 


2 


8.80 10"i 


2.33 10" 


1 


C4 


0.5 


6 


8.0 


1.128 


9.027 


[0.5, 1] 


47.0 




3.44 


2.78 10" 


3 


C4 


0.5 


6 


8.0 


1.128 


9.027 


[1,2] 


7.23 




37.37 


0.10 




C4 


0.5 


6 


8.0 


1.128 


9.027 


[2,3] 


6.96 10" 


2 


2.98 


1.22 10" 


1 


C4 


0.5 


6 


8.0 


1.128 


9.027 


[3, 14] 


2.30 10" 


2 


1.89 


2.33 10" 


1 


C5 


1.0 


6 


8.0 


2.257 


18.05 


[3, 14] 


2.30 10" 


2 


5.43 10-1 


4.15 10" 


3 


C6 


1.0 


4 


4.0 


1.263 


5.053 


[3, 14] 


2.30 10" 


2 


8.15 10-1 


1.45 10" 


1 
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Fig. 19. Angular distributions for C3 (see table [2)| for the mass bin 2 < Mjj < 3 TeV. Left: cross section in mb, right: cross 
section normalized to unit area 1 < x < 100. 
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Fig. 20. Angular distributions for C4 (see table [2| for the mass bin 1 < Mjj < 2 TeV. Left: cross section in mb, right: cross 
section normalized to unit area 1 < x < 100. 



is for each bin taken as ^/A^qcd,!:- Based on the fact that 
the theoretical uncertainty does not exceed 20% (see figure 
[iniright), and that experimental uncertainties reported by 
the Tevatron experiments are less than 11% [Ull], the sys- 
tematic error Ssys,i is taken to be 25% over the whole x 
range. Using (equation ((2T|) ') and the number of degrees 
of freedom, the probabihty P of having cJA^totai/dx/-^totai 
given the null hyptothesis is true, can be calculated. The 
null hypothesis of identity is rejected for P < 0.1. 

Let us focus on the mass bin 1 < Mjj < 2 TeV and 
consider only those worlds with n=6 extra dimensions and 
work with an integrated luminosity of 10 pb~^. We have 
used the recipe mentioned above to test several physics 
scenarios and the green coloured area in figure [23] shows 
for which model parameters — Mcs and Mg/Mp — the null 
hypothesis is rejected. This region is from now on called 
the region of discovery. 

As can be seen from the figure, large values of Mcff 
(Mcff > 1.5 TeV) and small values of Ms/Mp {Ms/Mp < 
1) fall outside the region of discovery. The reason is the 



absence of black hole creation because in that region the 
lower limit on the black hole mass is drastically increas- 
ing with decreasing Mg/Mp. See equation (PO)) and the 
discussion underneath. 

Above calculations have been done without parton show- 
ering, multiple interactions or hadronization. In figure [24l 
we compare the angular distributions with and without 
parton showers in the mass 3 TeV < Mjj bin for C3. The 
effect of parton showers is most visible at low x values. 



5 Conclusions 

We have discussed dijet angular distributions at ^/s = 14 
TeV. First we have performed a QCD study, then we have 
shown the distributions in a scenario with gravitational 
scattering and black hole formation in large extra dimen- 
sions. 

We have used two different programs, JETRAD and 
NLOJET++, for the calculation of QCD up to next-to- 
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leading order, and found that both programs generate sim- 
ilar results. The angular distributions at NLO are flatter 
than the Born calculations, especially at large values of x 
(X > 100), which is mainly caused by the fact that the 
running of as has less effect on NLO than on LO calcu- 
lations. Different jet algorithms tend to keep the shape of 
the distributions unchanged, but give a different normal- 
ization. 

We have investigated the systematic uncertainties com- 
ing from the choice of renormalization (fJ-B.) and factoriza- 
tion (fJ-p) scale and parton distribution function (PDF), 
and found that a change in fip and PDF mainly influences 
the normalization. On the other hand, a change in has 
an impact on both the normalization and the shape of the 
distributions; the distributions have a similar behavior at 
low X but tend to spread out as x increases. 

The effects on dijet angular distributions from gravi- 
tational scattering and black hole production, have been 
studied in the ADD scenario, with the extra requirement 
that the membrane on which the standard model fields 
are allowed to propagate, has a finite but small width. 
The model parameters are the fundamental Planck scale. 



the width of the membrane and the number of extra di- 
mensions, and it has turned out that the phenomenol- 
ogy is very much dependent on the fundamental Planck 
scale and the width of the membrane. For a fundamen- 
tal Planck scale of around 1 TeV and for a wide range 
of parameter settings, quantum gravity effects have been 
observed in mass bins above 1 TeV. For small widths of 
the membrane, gravitational scattering is the most im- 
portant process, while black hole formation dominates for 
wider membranes. In both cases, the effects mainly show 
up at small values of x- The same conclusions can be made 
for the normalized distributions. Using the shape of the 
distributions, rather than the absolute normalization, for 
X < 50, we have determined the region of parameter space 
that could be discovered with 10 pb~^ integrated lumi- 
nosity and a 25% systematic uncertainty in the mass bin 
1 < < 2 TeV. 

In conclusion, uncertainties from QCD that cannot be 
reduced by normalizing the distributions mainly show up 
at large values of x (x > 100), while effects from quantum 
gravity are mostly present at small values of x- Depending 
on the region of x under study, dijet angular distributions 
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Discovery potential for 1 < M;: < 2 TeV at 1 0pb'^ 
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Fig. 23. Discovery potential for the mass bin 1 < Aljj < 2 TeV at lOpb ^, assuming a 20% systematic uncertainty 



1e-06 



1e-07 r. 



1e-08 r 



1e-09 



— I 1 1 1 1 1 1 r 

QCD+GS+BH, no parton showers 

QCD+GS+BH, with parton showers 



J I I L 



10 20 30 40 50 60 70 80 90 100 
X 



0.01 r 



0.001 



— I 1 1 1 1 1 1 r 

QCD+GS+BH, no parton showers 

QCD+GS+BH, with parton showers 



J I I L. 



10 20 30 40 50 60 70 80 90 100 
X 



Fig. 24. Angular distributions in the mass bin 3 TeV < Mjj for C3, with and without parton showers. Left: cross section in 
mb, right: cross section normalized to unit area 1 < x < 100. 



can therefore either be used as a probe for new physics or 
as a test of QCD. 

The authors wish to thank Torbjorn Sjostrand and Leif Lonnblad 
for the many useful discussions. 
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